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Abstract. In many clusters of galaxies, the cooling time at the core of 
the intracluster medium is much less than the age of the system, suggest- 
ing that the the gas should continually lose energy by radiation. Simple 
thermodynamic arguments then require that the expected "cooling flow" 
should exhibit a specific spectroscopic signature, characterized by a dif- 
ferential emission measure distribution that is inversely proportional to 
the cooling function of the plasma. That prediction can be quantitatively 
tested for the first time by the Reflection Grating Spectrometer (RGS) ex- 
periment on XMM-Newton, which provides high resolution X-ray spectra, 
even for moderately extended sources like clusters. We present RGS data 
on 14 separate cooling flow clusters, sampling a very wide range in mass 
deposition rate. Surprisingly, in all cases we find a systematic deficit of 
low temperature emission relative to the predictions of the cooling flow 
models. However, we do see evidence for cooling flow gas at tempera- 
tures just below the cluster background temperature, Tq, roughly down 
to Tq/2. These results are difficult to reconcile with most of the possible 
explanations for the cooling fiow problem that have been proposed to 
date. We also present RGS data on the massive elliptical galaxy NGC 
4636. In this case, we detect evidence for resonance emission line scatter- 
ing of high oscillator strength Fe L-shell emission lines within the gaseous 
halo of the galaxy. The detection of that effect leads to very tight con- 
straints on physical conditions within the halo. However, here again, the 
expected signature of a cooling fiow is not detected, perhaps suggesting 
some fundamental uncertainty in our understanding of radiative cooling 
in low density cosmic plasmas. 
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1. Introduction 

The study of cooling flows in clusters of galaxies and the gaseous haloes of 
ellipticals has been a controversial area of research for nearly 30 years. As 
first emphasized independently by Cowie and Binney (1977), Fabian and Nulsen 
(1977), and Matthews and Bregman (1978), the coohng time for gas at the cores 
of these systems is often significantly less than the Hubble time. This implies 
that the gas should continually lose energy by radiation. A pressure driven 
cooling flow is thought to develop, causing the gas in the densest regions to 
slowly accrete onto the central dominant galaxies. The net result is a steady 
deposition of cold matter in the core of the cluster or massive galaxy. 

General support for this picture came from some of the first imaging and 
spectroscopic observations obtained with the Einstein Observatory. Images of 
many clusters indicate sharp peaks in the surface brightness distribution, as 
predicted by cooling fiow models. Where crude spectrophotometry has been 
available, one typically sees a marked softening of the spectrum in the central- 
most regions, again consistent with expectations. Spectroscopically, the presence 
of the cooling fiow should manifest itself as an excess of emission (over isother- 
mal models) at the lowest energies, or via the presence of soft X-ray emission 
lines inconsistent with the background gas temperature. The latter effect was 
detected with the Focal Plane Crystal Spectrometer on Einstein, for the gaseous 
halo around M 87 (Canizares et al., 1979, 1982). 

Nevertheless, subsequent observations led to refinements of the simplest 
cooling flow models, and also pointed to some inconsistencies. For example, 
for a homogeneous cooling flow, the mass deposition rate should be even more 
centrally peaked than is inferred from the X-ray brightness profiles (Johnstone et 
al. 1992). Typically, a mass deposition rate roughly proportional to the radius 
is derived. This was potentially explained by Nulsen (1986) as a consequence 
of thermal instability. Dense blobs form locally everywhere in the intracluster 
medium where the cooling time is less than the age of the system. The result is 
that the cooling gas resides in a multiphase medium, distributed over a rather 
large volume. However, this picture was criticized by Balbus and Soker (1989) 
and others, who argued that in the presence of the gravitational potential, the 
growth rates of the linear instability are too weak to allow the cooling blobs to 
form. 

A potentially more serious problem is that the end products of cooling are 
generally not seen in longer wavelength bands at the expected levels. For exam- 
ple, only a very small fraction of the cooling gas can form stars. H I absorption 
measurements have found no evidence for cold condensed clouds (O'Dca et al. 
1998), and while CO emission has been detected (Edge 2001), the amount of 
inferred molecular gas is a factor 10 below what is predicted by X-ray cooling 
estimates. 

Further, the first moderate resolution spectra obtained with the Einstein 
Solid State Spectrometer and with ASCA did not provide good quantitative 
fits to the simplest cooling fiow models. Although cool gas is clearly present 
in the spectrum, there is a dearth of emission at the lowest energies. This was 
interpreted as evidence for cold absorbing material intermixed with the X-ray 
emitting gas (White et al. 1991). However, if cold gas is indeed present in the 
intracluster medium, it remains unclear why it does not form stars. 
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Figure 1. The differential emission measure distribution of the iso- 
baric multi-phase cooling flow model. The peak of the emission of 
various Fe L ions are marked. By resolving emission lines from these 
Fe L ions with the RGS, the model can be tested in considerable detail. 



2. The Importance of High Resolution Soft X-ray Spectroscopy 

Independent of the details of the cooling flow process, the emitted X-ray spec- 
trum of the cooling gas can be robustly predicted using simple thermodynamic 
arguments: For an isobaric cooling flow, the differential luminosity emitted in a 
small temperature range dT is given by: 

dL = — kdT. 

2 fiuip 

dL can also be expressed in terms of the cooling function of the plasma: 

dL = A{T)dEM 

where dEM = n?dV is the differential emission measure. Thus: 

dEM _5 M k 
dT ~ 2 finip A{T) ' ^ ' 

The differential emission measure distribution predicted by this simple ar- 
gument is displayed in Figure 1. In calculating A(T), we have assumed one-third 
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cosmic abundances, which is typical of intracluster gas. Marked in the figure 
are the respective temperatures where each of the Fe L-shell ions becomes the 
dominant charge state. As can be seen, the Fe L complex very nicely samples 
the expected distribution. Hence by measuring the strengths of the various Fe 
L emission lines in a high resolution soft X-ray spectrum, we can accurately 
constrain cooling flow models. 

3. Relevance of XMM-Newton 

The first experiment capable of unambiguously making such measurements is 
the Reflection Grating Spectrometer (RGS) on the XMM-Newton Observatory. 
XMM-Ncwton incorporates three densely nested grazing incidence telescopes, 
collectively providing unparalleled area for imaging and spectroscopy. The RGS 
consists of arrays of grazing incidence reflection gratings, mounted immediately 
behind two of the three telescopes, which pick off roughly half the light in 
the beams, and disperse it to sets of dedicated charge-coupled device detec- 
tors (CCDs). For each of the two arrays, the gratings are mounted on an in- 
verted Rowland circle, which also includes the respective telescope focus and 
the readout strip of CCDs. The gratings are all identical, and they are mounted 
at the same graze angle with respect to the incident ray passing through each 
grating center. This configuration produces nearly stigmatic and aberration-free 
focussing at all wavelengths in the spectrum. For a more complete description of 
this experiment, see den Herder et al. (2002). The other experiments on XMM- 
Newton are the European Photon Imaging Camera (EPIC), arrays of CCDs 
mounted at the foci of the three telescopes (see Turner ct al. 2002 and Striider 
et al. 2002), and the Optical Monitor (OM), an optical/ultraviolet telescope 
co-aligned with with the X-ray instruments (Mason et al. 2002). 

The RGS bandpass (A = 5 to 38 A), was adjusted so as to optimally sample 
the broad range of Fe L-shell transitions, as well as the K-shell features associated 
with C, N, O, Ne, Mg, and Si. These are the most prominent emission lines 
expected for collisionally ionized gas, emitting in the temperature regime T = 
10^ to W^K, characteristic of the expected cooling flows in clusters and elliptical 
galaxies. The effective area over this band averages ~ 120 cm^, if we include 
both of the two identical instruments. For on-axis point sources, the spectral 
resolution is AA 0.06 A, varying only slowly as a function of wavelength. Like 
the transmission grating experiments on the Chandra Observatory, the RGS is 
a slitless spectrometer, which means that the resolution degrades for extended 
sources. Nevertheless, since the dispersion of the RGS gratings is very high, the 
degradation is fairly weak: 

AA pa O.lA X AO 

where A9 is the source size in arc-minutes. Thus, for only moderately extended 

sources like clusters and ellipticals, we still obtain fairly high resolution spectra. 
(In contrast, the spectral resolution of the Chandra experiments degrades by 
almost a factor of 100 for a source extent of 1 arc-minute.) 
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As a component of its guaranteed time observing program, the RGS instrument 
consortium observed a sample of 14 compact clusters and groups of galaxies, 
chosen to exploit the spectral sensitivity of the experiment (see Table 1). For a 
complete discussion of the analysis of these data sets, see Peterson et al. (2002). 
The sample includes clusters and groups at a range of temperatures, from 1 
to 10 keV. All of these systems were expected to host cooling-flows, with in- 
ferred mass deposition rates varying by three orders of magnitude (M ~ 1 to 
1000 Mq yr~^). Surprisingly, it was apparent from even the first of these obser- 
vations, that our data were incompatible with the robust spectral predictions of 
the standard cooling flow model. An example is illustrated in Figure 2, which 
displays the derived spectra for the massive cluster A 1835. The intracluster 
medium in this system has a background temperature kTo Ri 8 keV, and an 
inferred mass deposition rate M ~ 2300MQyr~-^ (see Peterson et al. 2001 for a 
more complete description of the analysis for this case.) In the top panel of the 
figure, we show the simulated cooling flow spectrum that we expected for the 
parameters appropriate to this source. The middle panel shows the comparison 
with the data. As can be seen, the highest ionization lines expected (Fe XXIII, 
Fe XXIV) are indeed detected in the data, however the model vastly overpredicts 
the intensity of the lower ionization Fe L-shell lines, as well as the O VIII and O 
VII K-shell transitions. The discrepancy indicates a dearth of cooler gas. In the 
bottom panel, we show the comparison to the same cooling flow model, but now 
with the predicted emission measure distribution truncated at kT = 2.7 keV, 
with no contributions from gas at lower temperatures. This somewhat arbitrary 
modification of the model does appear to provide a fairly good description of 
the data. 



Cluster 



Exposure (ks) Redshift 



Abeh 1835 
Abeh 665 
Abeh 1795 



36 0.2523 

20 0.1818 

40 0.0622 

38 0.0580 
26 0.0347 
54 0.0460 
29 0.0328 

39 0.0442 
33 0.0353 
36 0.0163 
50 0.0372 
42 0.0043 
48 0.0185 



Sersic 159-03 



2A0335+096 
Abcll 4059 
Abell 496 
MKW 3s 
Abcll 2052 
Abell 262 
Abeh 1837 
M87 



NGC 533 



Table 1. Clusters used in the RGS cooling-flow analysis. See Peter- 
son et al. (2002) for more details. 
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Figure 2. The failure of the coohng-flow model in the massive 
cooling-flow, Abell 1835. The top panel shows the spectral predic- 
tion for the cooling-flow model using the differential emission measure 
distribution in Figure 1 from the maximum temperature of 8 keV down 
to keV. The middle panel shows the comparison of the model (blue) 
with the data (red). The model clearly overpredicts a number of Fe 
L emission lines. The bottom panel shows the same model, but the 
emission below ~ 3 keV is ignored. Thus, there is evidence for cooling 
from 8 keV to 3 keV but no further. See Peterson et al. (2001) for 
more information. 
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Note that the data do not indicate the complete absence of a coohng flow. In 
particular, the detected Fe XXIII and XXIV lines are indicative of temperatures 
less than 4 keV, and cannot arise in the 8 keV gas that characterizes most of 
the intracluster medium. Interestingly, the intensities of these lines are roughly 
at the right level for the mass deposition rate inferred from the X-ray surface 
brightness. It is as if the cooling does occur as predicted, but the emergent 
luminosity somehow disappears from the soft X-ray band at temperatures below 
~ half the background temperature of the cluster. 

This same pattern is also found in all of the other clusters in our sample. 
A subset of the uicasured spectra are shown in Figure 3. Note that in all cases, 
Fe XXIII and Fe XXIV lines are detected prominently, while emission lines of 
lower ionization species (e.g. Fe XVII) are strongly suppressed. While O VIII 
Lyman-a (A = 18.97A) is always detected, its measured intensity is also always 
far weaker than predicted by the cooling flow models for reasonable values of 
the oxygen abundance (0/Fe > 0.2). 

5. Quantitative Spectral Fits 

To derive more quantitative constraints on the discrepancies with the cooling 
flow predictions, we have performed model fits to the data. Since clusters are 
extended sources, standard model fitting techniques cannot be easily employed 
with RGS data. Instead, we use a novel Monte Carlo technique developed ex- 
pressly for this application (Peterson, Jernigan, and Kahn 2002). The Monte 
Carlo correctly accounts for the off-axis behavior of the instrument response, 
and the effects of arbitrary selection cuts and transformations of the data. 

We generate photons within the cluster assuming a given spatial and spec- 
tral model. These are propagated through the instrument, yielding a distri- 
bution in the dispersion coordinate, the cross-dispersion coordinate, and the 
CCD pulse height. After applying the same cuts as are applied to the real 
data, the extracted spectrum from the Monte Carlo is compared to the data 
via a x^-statistic. An iterative technique is used to optimize the fit through 
the adjustment of the astrophysical parameters that specify the spatial/spectral 
model. 

We use a simple /3-model to characterize the distribution. The parameters 
are derived from the simultaneous EPIC images of the source. For the spectrum, 
we use the MEKAL (Mewe, Kaastra, and Liedahl 1995) model to generate the 
line and continuum emissivities. The shape of the spectrum is then determined 
by the differential emission measure distribution. 

Outside the "break radius" determined from the EPIC image, we assume 
that the gas is isothermal, and fit for the temperature and emission measure. 
Inside, we add in a cooling flow differential emission measure distribution, but 
we allow the mass deposition rate, M, to be normalized separately in different 
temperature bins: (^Tq to Tq), {jTo to ^Tq), (gTo to jTo), etc., where Tq is the 
background temperature. This choice of binning is arbitrary, but it works fairly 
well because the fractional ionization curves are roughly equally spaced in the 
logarithm of the temperature. The elemental abundances are also left as free 
parameters, but are assumed to be spatially uniform. The latter assumption 
does not have a big effect on the fits. 
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Figure 3. The failure of the coohng-flow model in a set of clusters. 
The blue points are the data, the green line is the cooling-flow model, 
and the red line is an empirical model where we have allowed the lu- 
minosity in various temperature bins to change. All clusters show the 
largest deviation with the cooling-flow model for Fe XVII (15 and 17 
A) and Fe XVIII (14.2 and 16 A). For further details, see Peterson et 
al. (2002). 
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Figure 4. Differential luminosity distribution of a sample of 14 clus- 
ters. In the left panel the differential luminosity is plotted against the 
absolute temperature, and in the right panel the differential luminos- 
ity is plotted against the fraction of the background temperature. The 
various clusters are plotted as different colors, and the points which are 
not upper limits are connected. Both plots are normalized to the pre- 
diction from the isobaric cooling-flow model, which would predict all 
points to track the line y=l. A number of detections and upper limits 
are well below that line, however, indicating the failure of the model 
at the lowest temperatures. The right plot reveals a more systematic 
trend in the failure of the model. See Peterson et al. (2002) for more 
details. 
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In general, the fits to the data with this model are quite good. This can 
be seen from Figure 3, where we show a representative sample of the data on 
individual clusters with the best fit models overlayed. The models appear to 
correctly account for the relative strengths of the Fe L lines, O VIII, Ne X, and 
Mg XII. 

In Figure 4, we plot the ratio of the derived differential luminosity for each 
of the temperature bins to the differential luminosity predicted by the simple 
cooling flow model. In Figure 4a, these points are plotted as a function of 
absolute temperature, and in Figure 4b, as a function of fractional temperature 
with respect to the background temperature for each individual cluster. Different 
clusters are represented by different color points, and the relative points for a few 
particular clusters are connected by lines for clarity. At the lowest temperatures, 
we typically derive upper limits, as indicated, but in all cases, we do derive clear 
detections of emission in at least two distinct temperature bins. 

The horizontal dashed lines in these figures represent what would be ex- 
pected if the cooling flow model were correct. While it is possible for the highest 
temperature points to lie above the line (due to contributions from the back- 
ground temperature gas), one would expect the detections and/or upper limits 
to be consistent with this line for all of the lower temperature points. As can 
be seen, in essentially all cases, at least some of the lower temperature points lie 
well below the line, indicating a dearth of cooler gas. When plotted as a function 
of absolute temperature (Figure 4a), there is no other obvious pattern to the 
data. However, when plotted as a function of fractional temperature (Figure 
4b), it is clear that all of the clusters behave in roughly the same way: there 
is a clear "drop-out" of emission at temperatures below ^Tq, where Tq is the 
background temperature. 

In fact, the various clusters appear to be consistent with a differential emis- 
sion measure distribution given by: 



with a ~ 1 to 2, as opposed to a = 0, as required by Equation (1). This 
empirical relationship is indicated by the dashed sloping lines in Figure 4b. 
Thus, in addition to the model being inconsistent with the data at the lowest 
temperatures, the entire shape of the differential emission measure distribution 
is inconsistent with the standard cooling-fiow model. 

6. Possible Explanations for the Deficit of Soft Emission 

Since the initial reports of some of the results we have presented here (Peterson 
et al. 2001, Kaastra et al. 2001, Tamura et al. 2001), a number of explanations 
have been proposed to account for the observed deficit of soft X-ray emission. 
Most of these work energetically, i.e. they incorporate additional heating or 
cooling mechanisms that can plausibly supply or remove energy at a rate com- 
parable to the missing soft X-ray luminosity. However, it still unclear whether 
any can successfully account for the systematic aspects of the observed discrep- 
ancies, e.g. the fact that the drop-out of emission always seems to occur at ^Tq 
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over a wide range in Tq. Below, we briefly discuss a few of the ideas which have 
been suggested: 

Heating. Heating by particle outflows associated with embedded active 
galactic nuclei (Rosner & Tucker 1989, Tabor Sz Binney 1993, David et al. 2001), 
or by subcluster mergers (Markevitch 2001), can supply sufficient energy deposi- 
tion to counteract the expected cooling. In the context of the multiphase model, 
the energy input must be distributed evenly everywhere over the thermally un- 
stable region. Otherwise local "pockets" of gas may still cool and condense. In 
addition, to suppress cooling in a temperature regime that is a fixed fraction of 
the background temperature, requires that the heating process be self-regulating, 
i.e. the heating rate has to nearly balance the cooling rate for every cluster, over 
a range of three orders of magnitude in cooling luminosity. This seems rather 
contrived, and, at the very least, such models require a fair degree of fine tuning 
to be considered viable. 

Non-standard cooling. A more efficient means of cooling than radiation in 
the soft X-ray regime could reduce the expected line emission. Simply varying 
elemental abundances docs not help much, since iron dominates both the cooling 
and the observed spectrum. Mixing the hot gas with cooler embedded gas to 
reradiate the energy at longer wavelengths (Begelman & Fabian 1990, Fabian 
et al. 2001) is a possibility, but it is by no means clear why this would always 
kick in at a fixed fraction of the background temperature. In addition, copious 
soft X-radiation can be produced by charge transfer at hot/cold interfaces, so 
it might not be that easy to truly suppress the soft X-ray emission line flux in 
such a scenario. 

Electron conduction. Heat conduction from hotter to colder phases could 
suppress further cooling as the dense condensations first begin to form (Tucker &i 
Rosner 1983, Stewart et al. 1984, Bertschinger Sz Mciksin 1986). The presence of 
tangled magnetic fields in the plasma can strongly suppress conduction down to 
levels where it becomes unimportant energetically (Chandran & Cowley 1998), 
but this issue continues to be debated theoretically (Narayan & Medvedev 2001). 
Observationally, electron conduction has been observed to be suppressed by 
factors ~ 100 or more at cold fronts in clusters (Markevitch 2000), so the viability 
of this explanation also remains suspect. 
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Figure 5. X-ray spectrum of the elliptical galaxy, NGC 4636. The 
spectrum is dominated by emission lines from Fe XVII, Fe XVIII, and 
O VIII. The red line is a spectral/spatial model using a Monte Carlo 
method to simulate the effects of resonant scattering. See Xu et al. 
(2002) for further details. 



7. RGS Spectroscopy of the Elliptical Galaxy NGC 4636 

NGC 4636 is an EO-1 giant elliptical located in the skirt region of the Virgo clus- 
ter. The gaseous halo associated with this galaxy has a total X-ray luminosity 
Lx ~ 5.7 X lO'^^ergs s~^, making it one of the three brightest ellipticals in the 
sky. The X-ray temperature inferred from ASCA measurements is kT ~ 0.76 
keV, with inferred elemental abundances ~ one third solar (Matsumoto et al. 
1997). At this temperature, and the implied electron density, the gas at the core 
of the halo should cool on timescales less than 10^ years, leading to a cooling 
flow with a mass deposition rate M 2 MQyi~^. 

The RGS spectrum of NGC 4636 is shown in Figure 5. As indicated, it 
is dominated by very strong K-shell line emission from N, O, Ne, and Mg, as 
well a rich Fe L-shell spectrum. In contrast to the cluster spectra, here we see 
very strong emission from the lower iron charge states, Fe XVII - XIX. This is 
a consequence of the much lower background temperature. The red line in the 
figure is our best fit thermal model for this spectrum, which involves a mean 
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temperature ^ 0.6 kcV and elemental abundances ranging from 0.5 to 1 in solar 
units. For further details on the spectral fit, see Xu et al. 2002. 

For this case, however, we can derive additional tight constraints on physical 
conditions in the halo core by examining the spatial profiles of the data in the 
cross-dispersion direction. In particular, the RGS provides a nearly stigmatic 
image of the source in the cross direction. The emission line images arc observed 
to be extended, and the particular profile that we observe for each line provides 
a clean indication of the distribution of the emission in that line within the halo. 

This is of interest for NGC 4636, because, for the inferred halo parameters, 
the medium should be partially optically thick to resonance scattering for the 
high oscillator strength Fe L transitions. We observe this effect directly in the 
data. In particular, the ratio of the Fe XVII 3s — 2p lines at 17.1 A to the Fe XVII 
3d — 2p lines at 15.0 A is found to be centrally peaked. Both of these sets of lines 
have similar cmissivity dependences on temperature, so this cannot be produced 
by a temperature gradient. It is due to the fact that the high oscillator strength 
15.0 A line scatters several times before exiting the halo, thereby broadening its 
spatial distribution. 

In fitting the data, we construct a Monte Carlo spatial/spectral model that 
is self-consistently constrained to fit both the spectrum and the cross-dispersion 
emission line profiles. The model is highly over-constrained. For example, the 
same density required to account for the line intensities, also yields the correct 
Fe L line optical depths to reproduce the observed spatial profiles. This can only 
work at one value of the distance, and the distance we infer is consistent with 
the known distance to the galaxy. 

So for this particular case, there is essentially no ambiguity about the phys- 
ical conditions in the emitting plasma. Nevertheless, the observed spectrum is 
still dramatically inconsistent with the cooling flow predictions. Specifically, at 
the implied mass deposition rate, the Fe XVII line intensities should be a fac- 
tor 10 brighter than observed! Similarly, the O VII line intensities are vastly 
over-predicted. 

Evidently, the multiphase cooling flow model also fails on the scale of an 
isolated elliptical galaxy, and one without a currently active nucleus. In our 
view, this single observation provides one of the most stringent challenges to our 
understanding of cooling in low density cosmic plasmas. 

8. Summary and Conclusions 

The flrst high spectral resolution X-ray observations of cooling flow clusters 
continue to challenge conventional models for these systems. Contrary to earlier 
beliefs, the deficit of soft emission is not due to excess embedded absorbing 
matter in the intracluster gas, but to a true lack of emission lines from the 
lowest temperature charge states. We find this problem in all the systems we 
have studied, from the most massive clusters to the gaseous haloes of individual 
elliptical galaxies. 

Empirically, our data appear consistent with the expression given by Equa- 
tion (2), where the emergent differential luminosity is proportional to the tem- 
perature raised to a power a ~ 1 to 2, rather than being independent of tem- 
perature, as expected from simple thermodynamic arguments. The solution to 
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this puzzle may represent some fundamental new physics in the cooling of low 
density, hot gas. If so, it may have profound implications for models of galaxy 
and structure formation in the early universe. 
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